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Abstract 

The annihilation of dark matter particles releases energy, ionizing some of the gas in the Universe. 
We investigate the effect of dark matter halos on reionization. We show that the effect depends on 
the assumed density profile, the particle mass, and the assumed minimum halo mass. For Navarro- 
Frenk- White halos and typical WIMPs, we find the effect to be quite small. However, light dark 
matter candidates in the MeV range can contribute significantly to reionization and can make an 
important contribution to the measured optical depth. This effect may be used to constrain light 
dark matter models. We also study the effect of varying the halo density profile on reionization. 
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I. INTRODUCTION. 



The absence of significant Lya absorption (the Gunn Peterson test[j|) in the spectrum 
of many quasars imphes that the Universe is highly ionized up to a redshift ~ 6. Observa- 
tions of the cosmic microwave background (CMB) by the Wilkinson Microwave Anisotropy 
Probe (WMAP)0 suggest that the Universe was reionized at a redshift 11, assuming full 
ionization at lower redshifts. Primordial stars and quasars are commonly believed to have 
played a dominant role in the reionization of the Universe. In this article, we investigate 
another possibility, namely whether radiation from the earliest dark matter halos could have 
contributed significantly to reionization. 

If dark matter particles in halos annihilate producing standard model particles, some of 
the released energy is absorbed by the gas, resulting in ionization. The effect of particle 
annihilation on the ionization of gas by a uniform distribution of dark matter was studied 
by [I, 0| • These authors however, concluded that WIMP dark matter is unlikely to have a 
significant effect on ionization. The effect of dark matter clumping was taken into account 
by , who modified the dark matter distribution by including a "boost factor" , and found 
that WIMP annihilation could be relevant to reionization, and could account for part of 
the observed optical depth, particularly for light dark matter candidates. In this article, we 
extend the analysis of |5|, by considering dark matter halos with a generalized Navarro- Frenk- 
White (NFW) profile. The number density of halos is determined by the Press-Schechter 
formula. We provide a detailed computation of the ionized gas fraction, and calculate the 
resulting optical depth, for different halo and particle parameters. We then show that our 
results may be used to constrain light dark matter candidates. 



A. Luminosity of halos. 

We fit each dark matter halo with a generalized NFW profile 



p(-) = . . .j:. , (1) 



(r/rs)° [1 + r/r,Y 

a = 1, j3 = 2 corresponds to the well known NFW profilej^. p{r) is the dark matter density 
at r, and ps and rg are constants. Let r2oo denote the radius at which the mean matter 
density p equals 200 times the cosmological average at the formation redshift zp, i.e., 

p{zf) = 200 p,n,,,{l + zp)' (2) 

where pc = 3Hq/SttG is the critical density, Hq is the Hubble parameter today and Q^n is the 
matter fraction, is the redshift at which a 1-cr fluctuation containing mass Mmin enters the 
non-linear regime[7|. The minimum mass Mmin is determined by the free streaming scale, 
which depends on the coupling of dark matter particles with standard model particles[3]. 
The mass in dark matter enclosed within r2oo 



Mdm(r200) = /dmM 
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■rloo /dm p{zf) (3) 



M = M(r2oo) is the halo mass. The concentration parameter C200 = fim/^s and /dm is the 
fraction of mass in dark matter which we set equal to fidm/^m = 0.8287|^. Note that we 
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have defined C200 at the time of halo formation zp. The luminosity of the halo per unit 
photon energy is then given by the expression: 

dL ^ , , , ^ c?AL 
— = Lo(M,C2oo) X -—^x 

(JjJb (XJj 

Lo = ^ fdrAnr'p'ir) (4) 

where x = E^/nidm and mdm is the particle mass, measured in units of energy, {cxav) is 
the annihilation cross section of the WIMPs times the relative velocity, averaged over the 
velocity distribution. We set {a^v) = 3 x 10~^^ cm^/s, in order to obtain fldmh'^ ~ 0.1 
today|8|. We will also assume that (cXaf) is independent of v. dN^/dE^ = m^^dN^/ dx 
is the number of photons released per annihilation per photon energy. Here, we restrict 
ourselves to tree level processes. We use the phenomenological result of for neutralinos, 
and express dN^/dx = ae~^^/x^'^, where a and b are constants for a particular annihilation 
channel. Averaging over the channels considered in [9j], we find a = 0.9, b = 9.56. The factor 
of 2 in the denominator accounts for the fact that two particles disappear per annihilation. 
Using Eq. [H we find 
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where e is a dimensionless cutoff scale, required to make the luminosity finite for a > 1.5. 
For the NFW profile, we have a = 1, /? = 2, e = 0, and 

^ 4oo [!-(! + C2oo)-^] 

3 [log(l + C200) - C200(l + C200)"^J 

II. IONIZATION OF GAS. 

The probability of ionization per unit volume at location s is given by 

p^^^ ^ [1 - Xionjs)] [1 + z{s)f a 

s is the distance travelled by a light ray from redshift infinity to redshift z (in the matter 
dominated era, z ^ 1) 

_ 2c 1 , . 

a is the scattering cross section of photons with gas atoms. Xion is the ionized fraction and 
rib [1 — Xion] is the comoving number density of bound atoms. For a mixture of 76% H and 
24% He, and assuming singly ionized He (the contribution of doubly ionized He is smaller, 
and is hence neglected here), we have 



0.82pc ^ 



b 



d% , < dL /r\\ 

r}p[s)n—- (9) 
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where 



0.76 



+ 



0.06 



0.82 13.6 eV 0.82 24.6 eV. 



0.07 eV" 



(10) 



i{s) is the number of ionizations per unit time per unit volume, at s. For simphcity, we have 
accounted for the different ionization potentials of H and He using /i defined in Eq. (TUl rather 
than by modifying the cross section. In any case, i{s) in Eq. Ois only an approximation, 
since we do not consider in detail, the scattering processes involved, t] is the fraction of 
the energy absorbed by the gas, resulting in ionization. We assume rj = 0.3, in accordance 
with [13]. We have assumed an ionization potential of 13.6 eV for Hydrogen and 24.6 eV for 
singly ionized Helium[ll|. f2b is the baryon fraction and nip is the proton mass. The cross 
section due to scattering of gas atoms by photons is the energy dependent Klein-Nishina 
cross section cr{E^) = f„{E^) 12(] (For photon wavelengths much smaller than the Bohr 
radius, we may neglect details of the bound structure while computing the scattering cross 
section. We therefore consider only electron-photon scattering.): 
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where y = E^/rrie, o"t is the Thomson cross section, and rrie is the electron mass. 

The fraction of the number of photons available at location s, having been emitted at s 



K[S, S 



exp 
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(12) 



In Eq. [121 "we have included both bound and ionized fractions in nb. Let dn/dM be the 
comoving number density of dark matter halos per unit halo mass, at a redshift z, given by 
the Press-Schechter formulahi 



dn 
dM 




2 Sc (1 + z) dah 



dM 



exp 



^g(l + ^) 
2al 



(13) 



Pm is the comoving matter density, is the variance of the density field smoothed over 
scale R 1^ (corresponding to a mass M) 



W{kR) 



dk k^P{k) 
~k 27r2 



\W{kR)f 



{kRf 



[sin kR — {kR) cos kR\ 



(14) 



P{k) is the matter power spectrum which takes the form given by Eisenstein and Hu[15 
normalized to = 0.8. Following |l6|, we set 6c = 1.28. We can now write an expression 
for the number of ionizations per unit volume per unit time, at redshift z: 



fz rl rip rao 

I{z)=Aixrj[l-x,„^]{l + zf -dz'{l + z')-^/^ dx^US 

Jzp Jo y/X Jm„ 



,M±-^U (15) 

In Eq. [m we have included all halos with masses above the cutoff mass M^^^. We then inte- 
grated over photon energies, and finally over redshift, from the formation redshift ZF(^min) 
to z. The factors of (1 + z) occur because we made a change from co-ordinate space to 
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redshift space, and also account for the redshifting of a photon as it propagates from z' to 
z. A is the dimensionless quantity 



and S'(z', z\ E^) is given by the expression 



A = 2=L ^ 3.5 X 10"=^ 



S = exp 



{(l + /)5/2_(l + ^)5/2| 



(16) 



(17) 



Let us now consider recombination into bound atoms. The recombination rate per unit 
volume 



0.76 0.06 



0.82 



0.82 



0.724 ^^3g-l 



S —1 

cm s 



(18) 



R{z) =nlxl^{l + z) 
ctH and ctHe are given by [13] 

an ~ 3.746 x 10-^^{T/eV) 
aue ~ 3.925 x IQ-^^iT/eVy^-^^^^ 

T is the gas temperature ^ 8 x 10^^[(1 + z)/61]^ eV. Using Eq. [15] and Eq. [181 we solve for 

•^ion (^) 



(19) 



I{z)-R{z) = n^il + z] 



.dx 



3 u-^ion 



dt 



where we have used the relation 

— Hq y f^ni dt 
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^ dz 

dz 



'l + z 
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1 + Z)J{1 + Zf + 



(20) 



(21) 



with the approximation Q\ <^ fim(l + z)^, valid at high redshifts. We neglect gas clumping 
in this work. The effect of gas clumping would be to increase the recombination rate. We 
however expect this effect to be small at redshifts z > 10. 



III. OPTICAL DEPTH. 



As mentioned earlier, most of the Hydrogen in the Universe is ionized up to redshift 
z = 6. We assume that Helium is singly ionized at 2; = 6 and doubly ionized 18| at 2: = 3. 
We then have the optical depth 



t{z < 6) = nb <7t 



0.88 
082 



fdz^{l + zf- fdz^{l + 
Jo dz J3 dz 



(22) 



which yields t{z < 6) = 0.04 where we have used Eq. [21] and the 5 year mean values of 
the WMAP experiment H] Hq = 71.9, = 0.742, J^^n = 0.258, = 0.04422. From the 
mean value of the WMAP measured optical depthf^] r = 0.087 ± 0.017, we note that an 
excess St ~ 0.047 needs to be provided by sources at redshift z > 6. We explore the dark 
matter models that can contribute this value 6t to the optical depth. In these models, 6t 
is provided solely by dark matter halos. 
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IV. RESULTS. 



We now present our results for Xion and 5t. We set {a^v) = 3 x 10^^^ cm^ s~^ and the 
density profile power law exponent [3 = 2 throughout. Our results thus depend upon a, 
the cutoff parameter e, the halo concentration parameter C200! the particle mass mdm, and 
the minimum halo mass Mmin- For simplicity, we treat C200 as a free parameter, and vary it 
independently of the halo mass. 



A. A generic WIMP. 

The earliest WIMPy halos are expected to form at a redshift z ~ 60, with masses M ~ 
[l9| . We now calculate the ionized fraction Xion using Eq. [201 and assuming a 
particle mass of mdm = 100 GeV, and Mmm = 1O~^M0. Xion depends on the power law 
exponent a, cutoff e, and concentration parameter C2oo- Fig[T]^a) shows function 
of redshift, for 4 different values of a, with C200 = 10, and e = 10~ . set equal 

to 1 for 2; > 6, in accordance with quasar studies. The NFW profile results in very little 
ionization, the peak ionization for z <Q being only 5 x 10~^. The ionized fraction increases 
with increase in a. For the case a = 1.9, we find a peak ionized fraction of 6%. Fig. Wljo) 
shows a^ion at ^ = 20 as a function of a, for different values of the cutoff parameter e, with 
C200 set equal to 10. Xion is almost independent of e for small a since the luminosity is finite 
as e — i> for a < 1.5. However, for larger values of a, Xion substantially increases as e ^ 0. 
For the extreme case of a = 2, e = 10~^, the ionized fraction is nearly 30%. Fig. [^c) shows 
the dependance of 6t on the halo concentration parameter 0200? for « = 1.6,1.8,2. The 
cutoff parameter e was set to 10~^. The solid red line indicates the maximum allowed value 



of 6t = 0.047. The solid black line shows the expected Planck sensitivity Tmin = 0.005[20 
As expected, the model with a = 2 produces the most optical depth and is excluded for 
C200 ~ 7. In all cases, the NFW model (a = 1) can be expected to be negligible in influencing 
reionization and the optical depth. 



B. Light dark matter candidates. 

We saw in the previous subsection that WIMPs with mass ~ 100 GeV do not cause sig- 
niflcant reionization, with the NFW density proflle. Let us now consider lighter dark matter 
candidates, with masses in the MeV range. The cosmology of MeV dark matter has been 
considered by many authors 21, [H, [H, 24]. Studies of the soft gamma ray background [23| 



and the 511 keV line|2J] can be used to constrain the particle mass. Here, we investigate 
an independent approach, namely using the WMAP constraint on r to place limits on the 
particle mass for a given cosmology. The form of dN^/dx depends on the particle physics of 
the theory. In this article, we assume the same dN^/dx as used previously for more massive 
WIMPs. Our results must therefore be treated with caution. 

We consider at flrst, halos with the NFW density proflle, i.e. with a = 1, e = 0. The 
results are sensitive to the chosen minimum halo mass Mmm- Afmin depends on the particle 
physics of the theory (in particular, the free streaming scale), and can be quite large (~ 



dwarf spheroidal galaxy) in certain models j22j. In those models, we do not expect to see a 



signiflcant effect on reionization. However, since Mmin is model dependent, we are justifled 
in treating it as a free variable, and we consider models in which Mmin varies from 10~® to 
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FIG. 1: (color online) Effect of neutralino dark matter halos on reionization and optical depth, (a) 
shows Xion for different values of the power law exponent a. zp ^ 57 for Mmin = 10~^Mq. Xion is 
set to 1 for z < 6. (b) shows the variation of Xion with the cutoff parameter e. The dependance on e 
is very slight for a < 1.5, but becomes important for larger values of a. (c) shows the contribution 
of dark matter halos to the optical depth 5t, for different values of C2oo- The solid red line indicates 
the maximum value of 6t consistent with the WMAP observation. The solid black line indicates 
the expected sensitivity of the Planck experiment. 
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IMq. We also assume that (o"af) = 3 x 10 cm^ s ^ independent of v 23 



Fig. [2]^a) shows the ionized fraction Xion, for mdm = 100, 10, and 1 MeV, for an assumed 
minimum mass Mmin = l^^o, and for C200 = 10. The halos are assumed to exist for 
z > zp, which in this case is zf{1 Mq) = 37.3. The ionization calculation was carried 
out to z = 6, beyond which, we set Xion = 1- We see that for MeV scale dark matter 
particles, the ionized fraction xion is significant (several percent) for NFW halos. There 
are two reasons that account for this increase. Firstly, a small particle mass results in a 
larger number density, and hence a higher probability of annihilation. Secondly, lighter dark 
matter particles produce photons with lower energy, resulting in a larger cross section for 
scattering. Fig. Mjo) shows the variation of xion with the minimum halo mass Mmin (note 
that zf depends on Mmin). Fig. [2]^c) shows the dependence on C2oo- 

Let us now consider the contribution of dark matter halos to the optical depth, for MeV 
dark matter. Fig. [3] shows 6t as a function of C200 for different values of mdm and Mmin- The 
four panels show the results for minimum halo masses Mmin = 1? 10~^, 10~^ and 10~^Mq. 
As expected, models with small Mmin produce the most 6t. The solid red line indicates 
the maximum allowed value of 6t = 0.047, while the solid black line indicates the expected 
sensitivity of the Planck experiment. The broken lines show 6t for different particle masses 
"^dm = 1, 10 and 100 MeV. We see that very light (~ 1 MeV) dark matter candidates result 
in too much St for models with Mmin < IO'^Mq, and are still disfavored for Mmin = 1 M©, 
except for very small values of C2oo- We are thus able to constrain the particle mass for given 
halos and particle physics models. 

Consider now, the effect of varying the density profile power law exponent a. Fig. H] 
shows 6t as a function of C200 for a particle mass m^m = 10 MeV, and M^in = 1 Mq. Three 
values of a are considered, namely a = 1, 1.3 and 1.6. The cutoff parameter e was set equal 
to 10""^. For 10 MeV candidates, we see that large a ~ 1.6 and small e ~ 10 ^ produce too 
much optical depth. As before, the solid red and black lines indicate the maximum 6t and 
the expected Planck sensitivity respectively. 



V. CONCLUSIONS. 

We fitted dark matter halos with a generic NFW profile parametrized by power law 
exponents a and (3, concentration parameter C200, and halo mass M. The luminosity of 
these halos was computed in terms of these parameters. {(Tav) = 3 x 10~^^ cm^s~^ was 
assumed to be independent of v. The Press-Schechter formalism was used to determine 
the number density of halos greater than a given minimum mass Mmin- We computed the 
ionization fraction as a function of redshift, produced by these dark matter halos. We then 
calculated the resulting optical depth and compared it to the WMAP measured value. 

We first analyzed the effect of 100 GeV WIMPs (e.g. the lightest supersymmetric particle 
or the lightest Kaluza-Klein particle) on reionization. We found that if the dark matter halos 
follow an NFW profile, the effect on reionization and the optical depth is negligible. To have 
a considerable effect, the density profile needs to be considerably more cuspy (a ^ 1.5), with 
small cutoff e ^ 10""^. 

We then looked at lighter dark matter candidates, in the 1 to 100 MeV range which 
were found to be more promising. We found that NFW halos have a significant effect on 
the ionized fraction, with particle masses mdm < 100 MeV. We analyzed the contribution 
of dark matter halos to the optical depth, for different values of the minimum halo mass 
Mmin- For small Mmin ~ 10^'^ Mq, we conclude that small particle masses mdm < 10 MeV 
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FIG. 2: (color online) Reionization due to halos composed of MeV scale dark matter particles, 
(a) shows Xion for different particle masses, zp was set to 37.3 for Mmin = ^Mq. xion = 1 for 
z < 6. Shown are results for particle masses mdm = 100, 10, and 1 MeV. (b) shows the variation 
of 

Xion with the minimum halo mass -Mrnin- (c) shows the dependence of Xion on the concentration 
parameter. 
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FIG. 3: (color online) Optical depth as a function of the concentration parameter, considering 
halos made of MeV dark matter. Each panel shows a different value of the minimum halo mass 
-^min- The solid red line shows the maximum allowed value of br. The solid black line indicates 
the expected sensitivity of the Planck experiment. The broken lines show the results for particle 
masses mdm = 100, 10, and 1 MeV. 

are not favored. For larger Mmin ~ 1 ""^dm ~ 10 MeV is allowed for small values of 
the concentration parameter C2oo- The value of Mmin is determined by the free streaming 
scale of the dark matter particles, which in turn depends on the scattering cross section 
of the dark matter particles with standard model particles in the early Universe. In order 
to determine Mmin, we would need to have a detailed understanding of the particle physics 
of the dark matter model. We can therefore use the predicted value of optical depth and 
future observations, to place constraints on MeV dark matter models. These limits provide 
information which may be used in conjunction with constraints on obtained by other 
means. 

Our model of reionization with MeV dark matter in NFW halos suggests a uniform, and 
gradual reionization scenario, rather than a "patchy" or "percolation" scenario wherein ion- 
ized bubbles form first, and then merge together. Future observations of the 21 cm emission 
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FIG. 4: (color online) Optical depth br as a function of the concentration parameter. Shown are 
cases a = 1, 1.3, 1.6. a = 1 corresponds to the NFW profile. The cutoff parameter e = 10~^, = 
10 MeV, and Mmin = IMq. The maximum allowed 5t = 0.047 is indicated by the solid red line. 
The expected Planck sensitivity is shown by the solid black line. 

of neutral H may help distinguish this model from the patchy reionization model. A better 
understanding of the reionization process could also be relevant to dark matter substructure 
theories. A significant ionized fraction of gas could have important consequences for the 
formation of molecular H2 which is the dominant coolant for primordial star formation. It 
would therefore be interesting to investigate whether early star formation could be influ- 
enced by a gradual reionization process such as the one we have described. Our results 
are also relevant to theories where early star formation may be affected by dark matter 
annihilation [25!]. Increased precision measurements of the CMB polarization would help us 
learn more about the reionization history. A detailed study of this epoch would lead to a 
better understanding of the first stars, and may also help us identify dark matter in the 
cosmic H-chamber. 
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